PARTICLE INTERFEROMETRY: 
NEW THEORETICAL RESULTS 



Ulrich Heinz 

Institut fur Theoretische Physik, Universitat Regensburg, 
D-93040 Regensburg, Germany 

By measuring hadronic single-particle spectra and two-particle correlations in 
heavy-ion collisions, the size and dynamical state of the collision fireball at freeze- 
out can be reconstructed. I discuss the relevant theoretical methods and their lim- 
itations. By applying the formalism to recent pion correlation data from Pb+Pb 
collisions at CERN we demonstrate that the collision zone has undergone strong 
transverse growth before freeze-out (by a factor 2-3 in each direction), and that it 
expands both longitudinally and transversally. From the thermal and flow energy 
density at freeze-out the energy density at the onset of transverse expansion can 
be estimated from conservation laws. It comfortably exceeds the critical value for 
the transition to color deconfined matter. 



1 Introduction 

In the last few years a large body of evidence has been accumulated that the 
hot and dense collision region in ultrarelativistic heavy ion collisions thermal- 
izes and shows collective dynamical behaviour. This evidence is based on a 
comprehensive analysis of the hadronic single particle spectra. It was shown 
that all available data on hadron production in heavy ion collisions at the 
AGS and the SPS can be understood within a simple model which assumes 
locally thermalized momentum distributions at freeze-out, superimposed by 
collective hydrodynamical expansion in both the longitudinal and transverse 
directionslircl The collective dynamical behaviour in the transverse direction is 
reflected by a characteristic dependence of the inverse slope parameters of the 
rn^-spectra ("effective temperatures") at small mi on the hadron massesGJ 
New data from the Au+Au and Pb+Pb systems □ support this picture and 
show that the transverse collective dynamics is much more strongly exhibited 
in larger collision systems than in the smaller ones from the first rounds of 
experiments. The amount of transverse flow also appears to increase mono- 
tonically with collision energy from GSI/SIS to AGS energies, but may show 
signs of saturation at the even higher SPS energy!! 

The extraction of flow velocities and thermal freeze-out temperatures from 
the measured single particle spectra relies heavily on model assumptionscl 
There have been alternative suggestions to explain the observed features of 
the hadron spectra without invoking hydro-nowBU The single-particle spectra 
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are ambiguous because they contain no direct information on the space-time 
structure and the space-momentum correlations induced by collective flow. In 
terms of the phase-space density at freeze-out ("emission function") S(x,p) 
the single-particle spectrum is given by EdN/d 3 p = f d 4 x S(x,p); the space- 
time information in S is completely washed out by integration. Thus, on the 
single-particle level, comprehensive model studies are required to show that a 
simple hydrodynamical model with only a few thermodynamic and collective 
parameters can fit all the data, and additional consistency checks are needed 
to show that the extracted fit parameter values lead to an internally consistent 
theoretical picture. The published literature abounds with examples demon- 
strating that without such consistency checks the theoretical ambiguity of the 
single particle spectra is nearly infinite. 

This is the point where Bosc-Einstein correlations between the momenta 
of identical particle pairs provide crucial new input. They give direct access to 
the space-time structure of the source and its collective dynamics. In spite of 
some remaining model dependence the set of possible model sources can thus be 
reduced dramatically. The two-particle correlation function C(q, K) is usually 
well approximated by a Gaussian in the relative momentum q whose width 
parameters are called "HBT (Hanbury Brown-Twiss) radii" . It was recently 
shown 11011 that these radius parameters measure certain combinations of the 
second central space-time moments of the source. In general thev mix the 
spatial and temporal structure of the source in a nontrivial wayji] and the 
remaining model dependence enters when trying to unfold these aspects. 

Collective dynamics of the source leads to a dependence of the 
radii on the pair momentum K; this has been known for m any ye> P 
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but was re cently quantitatively reanalyzed, both analytically LLH2H H I and 
numerically^ 3&3 The velocity gradients associated with collective expansion 
lead to a dynamical decoupling of different source regions in the correlation 
function, and the HBT radii measu re the size of the resulting "space-time re- 
gions of homogeneity" of the sourcell3liil around the point of maximum emis- 
sivity for particles with the measured momentum K . The velocity gradients 
are smeared out by a thermal smearing factor arising from the random mo- 
tion of the emitters around the fluid velocityu Due to the exponential decrease 
of the Maxwell-Boltzmann distribution, this smearing factor shrinks with in- 
creasing transverse momentum K± of the pair, which is the basic reason for 
the /^-dependence of the HBT radii. 

Unfortunately, other gradients in the source (for example spatial and tem- 
poral temperature gradients) can also generate a if-dependence of the HBT 
radiiSEl Furthermore, the pion spectra in particular are affected by resonance 
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decay contributions, but onl y a t s mall K±. This may also affect the HBT 
radii in a /-^-dependent wayj~~lcB The isolation of collective flow, in particu- 
lar transverse flow, from the if^-dependence of the HBT radii thus requires a 
careful study of these different effects. 

Our group studied this if-dependence of the HBT radii within a simple 
analytical model for a finite thermalized source which expands both longitudi- 
nally and transversally. For presentation I use the Yano-Koonin-Podgoretskii 
(YKP) parametrization of the correlator which, for sources with dominant 
longitudinal expansion, provides an optimal separation of the spatial and tem- 
poral aspects of the sourceElHil The YKP radius parameters are independent 
of the longitudinal velocity of the observer frame. Furthermore, in all thermal 
models without transverse collective flow, they show perfect Mj_-scaling (in 
the absence of resonance decay contributions). Only the transverse gradients 
induced by a non-zero transverse flow can break this Mj_-scaling, causing an 
explicit dependence on the particle rest mass. This allows for a rather model- 
independent identification of transverse flow from accurate measurements of 
the YKP correlation radii for pions and kaons. High-quality data should also 
allow to control the effects from resonance decays. 

Due to space limitations, I will be selective with equations, figures and 
references. A comprehensive and didactical discussion of the formalism and 
a m o re o vtnn si ve selection of numerical examples can be found in the lecture 
notesO to which I refer the reader for more details. 



2 Spectra and emission function 

2.1 Single-particle spectra and two-particle correlations 

The covariant single- and two-particle distributions are defined by 



dN 
d 3 p 

dN 



p iiP) = E ~ =E(a+a p ), (1) 



Pi (Pa ,Pb) = E a E b = E a E b (a+ a o+ a Pb d Pa ) , (2) 

where dp (d p ) creates (destroys) a particle with momentum p. The angular 
brackets denote an ensemble average, 

(6)=tr(pO), (3) 

where p is the density operator associated with the ensemble. (When talking 
about an ensemble we may think of either a single large, thermalized source, 
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or a large number of similar, but not necessarily thcrmalizcd collision events.) 
The single-particle spectrum is normalized to the average number of particles, 
(N), per collision, 



(4) 



while the two-particle distribution is normalized to the number of particles in 
pairs, (N(N — 1)), per event: 



/ 



^^P 2 (Pa, Pb ) = (N(N-l)). (5) 
The two-particle correlation function is defined as 

If the two particles are omi tted independently and final state interactions are 
1 7 

neglected one can proveLLil a generalized Wick theorem 

fir,,) I 2 

C(p a ,p b ) = l± p 4'. (7) 

{ap a a P J {ap b a Pb ) 

Note that the second term is positive definite, i.e. the correlation function 
cannot, for example, oscillate around unity. This is no longer true if final state 
interactions are included (see below). 

From now on I will assume that the emitted particles are bosons, and for 
convenience I will call them pions, although nearly everything below applies 
equally well to other bosonic particles. 



2.2 Source Wigner function and spectra 

In the language of the covariant current formalismlll the source of the emitted 
pions can be described in terms of classical currents J(x) which act as classical 
sources of freely propagating pions. These currents represent a parametrization 
of the last collision from which the free outgoing pion emerges. Very helpful 
for the following will be the so-called "emission function" S(x, K)ur^i 

It is the Wigner transform of the density matrix associated with the classical 
source amplitudes J(x). This Wigner density is a quantum mechanical object 
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defined in phase-space (x,K); in general it is real but not positive definite. 
But, when integrated over x or K it yields the classical (positive definite and 
real) source density in momentum or coordinate space, respectively, in exactly 
the same way as a classical phase-space density would behave. Furthermore, 
textbooks on Wigner functions show that their non-positivity is a genuine 
quantum effect resulting from the uncertainty relation and are concentrated at 
short phase-space distances; when the Wigner function is averaged over phase- 
space volumes which are large compared to the volume (27r?i) 3 of an elementary 
phase-space cell, the result is real and positive definite and behaves exactly like 
a classical phase-space density. 

The emission function S(x, K) is thus the quantum mechanical analogue 
of the classical phase-space distribution which gives the probability of finding 
at point x a source which emits free pions with momentum K. It allows to 
express the single-particle spect ra an d t wo-particle correlation function via the 
following fundamental relationsElll^l'^ 

E P1T- = fd 4 xS(x,p), (9) 



d 3 p 

\I d 4 xS(x,K)e l i- x \ 2 
C(q,K) - 1 + -j — — - - j - j — — R _i q y (10) 

For the single-particle spectrum (||), the Wigner function S(x,p) on the r.h.s. 
must be evaluated on-shell, i.e. at p° = E p = y/ m 2 + p 2 . For the correlator 
( |l0| ) we have defined the relative momentum q = p a — pb, q° = E a — Eb 
between the two particles in the pair, and the total momentum of the pair 
K = {p a +Pb)/2, K° — (E a + Eb)/2. Of course, since the 4-momenta p a .b of 
the two measured particles are on-shell, p® — Ei — y/ m 2 + pf , the 4-momenta 
q and K are in general off-shell. They satisfy the orthogonality relation 

q-K = 0. (11) 

Thus, the Wigner function on the r.h.s. of Eq. ( |To| ) is not evaluated at the 
on-shell point K° = Ek- This implies that for the correlator, in principle, we 
need to know the off-shell behaviour of the emission function, i.e. the quantum 
mechanical structure of the source. Obviously, this makes the problem appear 
rather untractable! 

Fortunately, nature is nice to us: the interesting behaviour of the correlator 
(its deviation from unity) is concentrated at small values of |q|. Expanding 
K° = (E a + E b )/2 for small q one finds 

A - o =M i+ 8i +o G?)H' f - (12) 
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Since the relevant range of q is given by the inverse size of the source (more 
properly: the inverse size of the regions of homogeneity in the source - see 
below) , the validity of this approximation is ensured in practice as long as the 
Compton wavelength of the particles is small compared to this "source size" . 
For the case of pion, kaon, or proton interferometry for heavy-ion collisions 
this is true automatically due to the rest mass of the particles: even for pions 
at rest, the Compton wavelength of 1.4 fm is comfortably smaller than any 
typical nuclear source size. This is of enormous practical importance because it 
allows you essentially to replace the source Wigner density by a classical phase- 
space distribution function for on-shell particles. This provides a necessary 
theoretical foundation for the calculation of HBT correlations from classical 
hydrodynamic or kinetic (e.g. cascade or molecular dynamics) simulations of 
the collision. 

If the single-particle^ spectrum is an exponential function of the energy 
then it is easy to prove 111 that one can replace the product of single-particle 
distributions in the denominator of ( [To| ) by the square of the single-particle 
spectrum evaluated at the average momentum K: 



J cftxe 1 ** S(x,K) 



J d 4 xS(x,K) 



= 1 + | (e lq - x ) | . (13) 

The deviations from this approximation are proportional toJ;he curvature of 
the single-particle distribution in logarithmic representation!!! They are small 
in practice because the measured single-particle spectra are usually more or 
less exponential. In the second equality of (|l^) we defined (. . .) as the average 
taken with the emission function; due to the -fT-dependence of S(x, K) this 
average is a function of K. This notation will be used extensively in the 
following. 

The fundamental relations (|J) and (jl^) rcsp. ( jl3| ) show that both the 
single-particle spectrum and the two-particle correlation function can be ex- 
pressed as simple integrals over the emission function. The emission function 
thus is the crucial ingredient in the theory of HBT interferometry: if it is 
known, the calculation of one- and two-particle spectra is straightforward (even 
if the evaluation of the integrals may in some cases be technically involved) ; 
more interestingly, measurements of the one- and two-particle spectra provide 
access to the emission function and thus to the space-time structure of the 
source. This latter aspect is, of course, the motivation for exploiting HBT in 
practice. In my talk I will concentrate on the question to what extent this ac- 
cess to the space-time structure from only momentum-space data really works, 
whether it is complete, and (since we will find it is not and HBT analyses 
will thus be necessarily model-dependent) what can be reliably said about the 
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extension and dynamical space-time structure of the source anyhow, based on 
a minimal set of intuitive and highly suggestive model assumptions. 



2.3 Final state interactions (FSI) 

Equation (|l^) reflects the absence of final state interactions (free propagation 
after emission) by the appearance of the plane wave factor e lq ' x under the 
integral of the exchange term in the two-particle cross section: 



P2(Pa,Pb) = J d 4 Xl d 4 X 2 [s(x U K+!Z) S(x 2l K- 



±e ig-(*i-»a) S(x u K)S(x 2 ,K) 



(14) 



In practice particle interferometry is done with charged particle pairs which 
suffer a long-range Coulomb final state repulsion on their way out to the detec- 
tor. In addition, there may be strong final state interactions, e.g. in proton- 
proton interferometry where there is a strong s-wave resonance just above the 
two-particle threshold. In this case Eq. (14) must be replaced bycil 



P2(Pa,Pb) = j ' d A xd A yS{x + ^, Pa ) S(x-%,p b ) 



0(y°) |<^/ 2 (y-W)| 2 + 0(-y°) \<l> q /2{y-v a y )\ 



± 



d 4 x d A y S (x 



IK 



x &- q/2 (y-vy a ) <t> q /2{y-vy°). 



(15) 



Here v=K / Ek , v a =p a / Ek, Vb=Pb/ Ek are (to quadratic accuracy in q) the 
velocities of the particles with momentum p a , pb, respectively, and 4> q /2(f) 
is an FSI distorted wave with asymptotic relative momentum q/2, evaluated 
at the two-particle relative distance r. Upon replacing the latter by plane 
waves (15) turns into (E3). The FSI distorted waves can be calculated by 
solving a non-relativistic Schrodingcr equation for the relative motion which 
includes the FSI potential in the rest system of the pair (where K = v = 0) . 
Eq. ( |l5| ) represents a non-relativistic Galilei-transformation of the result from 
the pair rest frame to the frame in which p a and p^ are measured; therefore it 
is only valid in observer frames in which the pair moves non-relativistically. In 
order to evaluate Eq. fllq) one must therefore first transform the 4-momenta 
p a ,b to such a frame (best directly into the pair rest frame). The momentum 
argument q of the FSI distorted waves <j> is then the difference between the two 
spatial momenta in that frame, and their space-time argument y — Viy° is the 
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relative distance of the two particles in that frame at the time when the second 
particle is emitteacA Since the latter depends not only on the time difference y° 
between emission points, but also on the velocity of the first emitted particle, 
these arguments depend on the momentum argument of the emission function 
associated with the first emitted particle. The two terms ~ 6(±y°) in the direct 
term reflect the two possible time orderings between the emission points. 



2.4 Implementation in event generators 

Equations (|TJ) and (|l5|) can be implemented into event generators, following 
the procedure given rnEiMl 

For the direct term one selects all pairs (i,j) with pi — p a , pj — p within 
a given numerical accuracy (bin width) which is essentially dictated by event 
statistics. Each pair is multiplied with a weight given by the corresponding 
probability density \<j> q /2\ 2 of the FSI distorted wave. The latter must be 
evaluated in a frame in which the pair moves non-relativistically, best in the 
pair rest frame where K — (p a +p ) /2 = 0. (Then Ek = m, and the velocities 
v,v a ,b reduce to their usual nonrelativistic definition.) From the space-time 
coordinates Xi , Xj of the pair in the event generator frame one calculates the 
distance y*j between the two particles in the pair rest frame at the time when 



the second particle is produced. One then computes \<j> q * /2(y*j)\ an d weights 
the selected pair with this number. In this expression q* is the spatial 
relative momentum between the two particles in the pair rest frame which must 
be computed from p a ,Pb in the event generator frame. In the absence of FSI, 
the corresponding weight is simply 1. The complete direct term is obtained by 
summing over all such pairs. 

For the exchange term, the selection of pairs and weights is a little less 
intuitive!!! One selects all pairs with pi = pj — K (i.e. on-shell particles 
(!) with Pi = pj — K and = Ej — Ek), again within the same numerical 
accuracy (bin width) as above. From the production coordinates Xi , Xj one 
again computes the spatial distance y*j between the two particles in the pair 
at the time of emission of the second one, in the pair rest frame K = 0. This 
distance is used to compute the weight </>*L q ,/ 2 (y*j)(l> q ,/ 2 (yij) for this pair. The 
value of q* here is the same as above in the direct term, i.e. it is computed 
from p a and pb by transforming into the pair rest frame, not from pi = pj = K. 
Without FSI, the corresponding weight ?1 is cos(q* • y*j)- The full exchange 
term is obtained by summing over all such pairs. Note that this selection 
of pai rs a nd weights differs from previously applied algorithms which were 
shown to yield wrong results for sources with very strong x-p correlations. 

Finally one must normalize the correlator by the product of single particle 
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spectra, 

Pi(p a )Pi(p b )= I d 4 xS(x, Pa ) [ d 4 yS(y, Pb ) 



(16) 



This normalization is best obtained from the pairs selected for the direct term 
above by multiplying them with unit weights. 

One may object to the use of event generators for the emission function 
because they fix particle momenta and coordinates simultaneously and thus 
violate the uncertainty principle. One can generate from an event genera- 
tor a quantum mechanically consistent Wigner density S(x,p) by folding the 
event generator output with minimum uncertainty wave packetsEHI The corre- 
sponding quantum mechanically consistent algorithm^ for computing single- 
and two-particle spectra is easily generalized to include FSI effects, by sim- 
ply replacing the factors 1 and cos(q* • y*j) in the direct and exchange terms, 
respectively, by the correct FSI weights as discussed above. 

At this point I will drop the discussion of final state interactions; the rest 
of the lecture will deal only with the case of free particles, assuming (carelessly) 
that appropriate Coulomb corrections of the data have already been done by 
the experimentalists. 



2.5 The mass-shell constraint 

Expressions (10 13) show that the correlation function is related to the emis- 
sion function by a Fourier transformation. At first sight this might suggest that 
one should easily be able to reconstruct the emission function from the mea- 
sured correlation function by inverse Fourier transformation, the single particle 
spectrum (^|) providing the normalization. This is, however, not correct. The 
reason is that, since the correlation function is constructed from the on-shell 
momenta of the measured particle pairs, not all four components of the relative 
momentum q occurring on the r.h.s. of ( |l3|) are independent. They are related 
by the "mass-shell constraint" ( |ll|) which can, for instance, be solved for q°: 

q °=f3-q with /3 = -^«|-. (17) 

f3 is (approximately) the velocity of the cm. of the particle pair. The Fourier 
transform in (|l^) is therefore not invertible, and the reconstruction of the 
space-time structure of the source from HBT measurements will thus always 
require additional model assumptions. 
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It is instructive to insert ( |l7| ) into (|l3|): 



C(q,fST)«l + 



/ d A x exp(iq-(x-pt)) S(x,K) 



jd A xS(x,K) 



(18) 



This shows that the correlator C(q, K) actually mixes the spatial and tempo- 
ral information on the source in a non-trivial way which depends on the pair 
velocity f3. Only for time-independent sources things seem to be simple: the 
correlator then just measures the Fourier transform of the spatial source distri- 
bution. Closer inspection shows, however, that it does so only in the directions 
perpendicular to (3 since the time integration leads to a ^-function 5{f3-q): 



lim 



J_ T dt exp (— i q f3 t) 



J_ T dt 



r lim^%-/3) 



(19) 



This implies that there are no correlations in the direction parallel to the pair 
velocity (3 (which will be called the "outward" direction below), i.e. C — 1 for 
<Zout 7^ 0. The width of the correlator in this direction vanishes! This should 
puzzle you: wouldn't you have thought that the width of the correlator in 
the "outward" direction is inversely related to the source size in that direction 
(which is, of course, perfectly finite)? As we will see in the next subsection this 
unexpected behaviour is just another consequence of the mixing of the spatial 
and temporal structure of the source in the correlator: The width parameter 
of the correlator in the "outward" direction receives also a contribution from 
the lifetime of the source which in this case diverges, leading to the vanishing 
width of the correlator. 

It is instructive to look at the problem also in the following way: If one 
rewrites Eq. ([l8|) in the pair rest frame where K = and hence q a = 0, one 
obtains 

C(q,K)-l = J d 3 r cos(q-r)S Te i(r;K) (20) 



where 



with 



S iel (r;K) = J d 3 Rs K (R+tr) s K (R-\r), 



(21) 



(22) 



/ d 4 x'S(x',K) ' 

is the time-integrated normalized relative distance distribution in the source. 
The latter can, in principle, be uniquely reconstructed from the measured cor- 
relator by inverting the cosine- Fourier transform (poj). But since it gives 
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only the time integral of the relative distance distribution for fixed pair mo- 
mentum K in the pair rest frame, no direct information on the time structure 
of the source is obtainable! Only by looking at the result as a function of 
K, which, as I will show, brings out the collective dynamical features of the 
source, can one hope to unfold the time-dependence of the emission function. 
It is clear that this will be only possible within the context of specific source 
parametrizations. 

2.6 K- dependence of the correlator 

We have seen that in general the correlator is a function of both q and K. Only 
if the emission function factorizes in x and K, S(x,K) = F(x) G(K), which 
means that every point x in the source emits particles with the same mo- 
mentum spectrum G(K) (no "x-if-correlations" Vthe if-dependence in G{K) 
cancels between numerator and denominator of (|13| ) , and the correlator seems 
to be iC-independent. However, not even this is really true: even after the 
cancellation of the explicit independence G(K), there remains an implicit in- 
dependence via the pair velocity f3 f=a K / Ek in the exponent on the r.h.s. of 
Eq. (|l^)! Only if both conditions, factorization of the emission function in x 
and K and time-independence of the source, apply simultaneously, the correla- 
tion function is truely ^-independent (because then the /3-dependence resides 
only in the 5- function ([l9|)). 

The only practical situation which I know where this occurs and a K- 
independent correlation function should thus be expected is in HBT interfer- 
ometry of stars for which the method was invented Ej It is hard to believe that 
this complication in the application of the original HBT idea to high-energy 
collisions went nearly unnoticed for more than 20 years and was stumbled 
upon more or less empirically by Scott Pratt in his pioneering work on HBT 
interferometry for heavy-ion collisions □ only in 1984! 

If one parametrises it by a Gaussian in q (see below) this means that 
in general the parameters ("HBT radii") depend on K. Typical sources of 
x-K correlations in the emission function are a collective expansion of the 
emitter and/or temperature gradients in the particle source: in both cases 
the momentum spectrum ~ exp[— p-u{x) /T(x)] of the emitted particles (where 
u^lx) is the 4- velocity of the expansion flow) depends on the emission point. In 
the case of collective expansion, the spectra from different emission points are 
Doppler shifted relative to each other. If there are temperature gradients, e.g. 
a high temperature in the center and cooler matter at the edges, the source 
will look smaller for high-momentum particles (which come mostly from the 
hot center) than for low-momentum ones (which receive larger contributions 
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also from the cooler outward regions). 

We thus see that collective expansion of the source induces a if-dependcncc 
of the correlation function. But so do temperature gradients. The crucial 
question is: does a careful measurement of the correlation function, in par- 
ticular of its X-dependence, permit a separation of such effects, i.e. can the 
collective dynamics of the source be quantitatively determined through HBT 
experiments? We will see that this is not an easy task; however, with suffi- 
ciently good data, it should be possible. In any case, the if-dependence of the 
correlator is a decisive feature which puts the HBT game into a completely 
new ball park. Two-particle correlation measurements which are not able to 
resolve the if-dependence of the HBT parameters arc, in high energy nuclear 
and particle physics, of very limited use only. [Unfortunately, this applies to all 
published HBT data from pp and e + e~ collisions. In my opinion, a renewed 
investigation of two-particle correlations from pp and e + e~ collisions, using 
the powerful new tool of multidimensional, differential HBT analysis, should 
be a high priority project - as it is, we have practically nothing with which to 
compare our heavy-ion results in a meaningful way.] 

3 Model-independent discussion of HBT correlation functions 

3.1 The Gaussian approximation 

The most interesting feature of the two-particle correlation function is its half- 
width. Actually, since the relative momentum q = pi — p 2 has three Cartesian 
components, the fall-off of the correlator for increasing q is not described by 
a single half- width, but rather by a (symmetric) 3x3 tensorl! which describes 
the curvature of the correlation function near q = 0. We will see that in 
fact nearly all relevant information that can be extracted from the correlation 
function resides in the 6 independent components of this tensor. This in turn 
implies that in order to compute the correlation function C it is sufficient to 
approximate the source function S by a Gaussian in x which contains only 
information on its space-time moments up to second order. 

Let us write the arbitrary emission function S(x, K) in the following form: 

S(x, K) = N(K) S(x(K),K) e-\^ [K) B ^ [K) *" {K) + SS(x, K) , (23) 

where we adjust the parameters N(K), x^(K), and B^K) of the Gaussian 
first term in such a way that the correction term 6S has vanishing zeroth, first 
and second order space-time moments: 

J d^xSSix.K) = J d A xx»5S{x,K) = J d 4 x x fl x 1/ SS(x, K) = . (24) 
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This is achieved by choosing 

, v dN det 1/2 B U JK) , . 

x»(K) = (a"), (26) 
(B- l )^{K) = (x^)^((x-x)^(x~x) u ). (27) 

The (if-dependent) average over the source function (. . .) has been defined 
in Eq. (|l3|). The normalization factor ( p5f ) ensures that the Gaussian term in 
( p3| ) gives the correct single-particle spectrum (^) ; it fixes the normalization on- 
shell, i.e. for K a = Ek, but as we discussed this is where we need the emission 
function also for the computation of the correlator. x(K) in ( p6| ) is the centre 
of the emission function S(x, K) and approximately equal to its "saddle point" , 
i.e. the point of highest emissivity for particles with momentum K. The second 
equality in ( p7j) defines x as the space-time coordinate relative to the centre of 
the emission function; only this quantity enters the further discussion, since, 
due to the invariance of the momentum spectra under arbitrary translations of 
the source in coordinate space, the absolute position of the emission point is 
not measurable in experiments which determine only particle momenta. Since 
x{K) is not measurable, neither is the normalization iV(.fr)liil as its definition 
( p5| ) involves the emission function at x{K). Finally, Eq. (|27j ) ensures that 
the Gaussian first term in ( |23| ) correctly reproduces the second central space- 
time moments (x^Xv) of the original emission function, in particular its r.m.s. 
widths in the various space-time directions. 

Inserting the decomposition ( |23| ) into Eq. ( [l3| ) we obtain for the correlation 
function 

C(q, K) = 1 + exp [-g" q v (K)] +SC{q,K). (28) 

The Gaussian in q results from the Fourier transform of the Gaussian contribu- 
tion in (HH); the last term SC receives contributions from the second term SS 
in ( ^ ) which contains information on the third and higher order space-time 
moments of the emission function, like sharp edges, wiggles, secondary peaks, 
or non-Gaussian tails in the source. It is at least of fourth order in q, i.e. the 
second derivative of the full correlator at q — is given exactly by the Gaussian 
in ( ^8| ) . Please note that the exponent of the correlator contains no term linear 
in q; since the correlator must be symmetric under q — > —q because it does 
not matter which of the two particles of the pair receives the label 1 or 2, a 
linear g-dependence could only arise in the form exp(— i?|q|). The only type 
of emission function yielding such a g-dependence of the correlator would be a 
spherically symmetric Lorentzian. Any emission function which at large x falls 
off faster than l/x 2 results in the leading Gaussian behaviour (EH) instead. 
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This settles, in my opinion, the old issue whether Gaussian or exponential fits 
of the correlation function should be preferred. 

[In the past it has repeatedly been observed that the correlation data 
appear to be better fit by exponentials than by Gaussians. However, as far as 
I know, this happened always when one tried to fit the correlator as a function 
of the single Lorentz invariant variable Qf nv = (q ) 2 — q 2 ■ Contemplating the 
structure of Eq. ( p8[ ) one realizes that such a fit does not make sense: the 
generic structure of the exponent, —q IJ 'q u {x l j l x v ), tells us that the term (q a ) 2 
should come with the time variance of the source while the spatial components 
(q 1 ) 2 should come with the spatial variances of the source. Since all variances 
are positive semidefinite by definition, it does not make sense to parametrize 
the correlation function by a variable in which (q ) 2 and q 2 appear with the 
opposite sign! Such a fit could only work if the time variance and all mixed 
variances would vanish identically, and all three spatial variances were equal. 
This is certainly not the general case in nature. The good exponential fits of 
the correlation functions from pp and e + e~ collisions are thus, in my mind, 
purely accidental and an empirical curiosity without physical meaning. The 
variable Qi nv should not be used for fitting HBT data.] 

Please note also that Eq. ( p8| ) has no factor i in the exponent. If the 
measured correlator is fitted by a Gaussian as defined in (|2^), its q- width can 
be directly interpreted in terms of the r.m.s. widths of the source in coordinate 
space. Any remaining factors of \pl, \f?>, or \fh (which you can sometimes find 
in the literature) are due to reexpressing the r.m.s. width of the source in 
terms of certain other width parameters chosen for the paramctrization of the 
source in coordinate space. The confusion connected with such factors is easily 
avoided by always expressing the source parametrization directly in terms of 
r.m.s. widths. 

Eqs. ( |23| ) and ( p8| ) would, of course, not be useful if the contributions from 
5S and 8C were not somehow small enough to be neglected. This requires a 
numerical investigation. It was shown 113 that in typical (and even in some 
not so typical) situations SS has a negligible influence on the half width of 
the correlation function. It contributes only weak, essentially unmeasurable 
structures in C(q, K) at large values of q. The reader can easily verify this 
analytically for an emission function with a sharp box profile; the resultsil for 
the exact correlator and the one resulting from the Gaussian approximation 
( p3| ) differ by less than 5% in the half widths; the exact correlator has, as a 
function of q, secondary maxima with an amplitude below 5% of the value 
of the correlator at q = 0. We have checked that similar statements remain 
even true for a source with a doughnut structure, i.e. with a hole in the 
middle, which was obtained by rotating the superposition of two 1-dimensional 
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Gaussians separated by twice their r.m.s. widths around their center. The 
only situation where these statements require qualification is if the correlator 
receives contributions from the decay of long-lived resonances; unfortunately, 
this is of relevance for pion interferometry as will be discussed in Sec. refsec6. 

From Eq. (|2^) we conclude that the two-particle correlation function mea- 
sures the second central space-time moments of the emission function. That's 
it - finer features of its space-time structure (edges, wiggles, holes) cannot 
be measured with two-particle correlations, but require the analysis of three-, 
four-, . . . , many-particle correlations. 1^3 The variances (x^Xu) are in general 
not identical with our naive intuitive notion of the "source radius" : unless the 
source is stationary and has no x- K -correlations at all, the variances depend 
on the momentum K of the pair and cannot be interpreted in texms of sim- 
ple overall source geometry. Their correct physical interpret at ionB'B^Hjj is in 
terms of "lengths of homogeneity" which give, for each pair momentum K , the 
size of the region around the point of maximal emissivity x(K) over which the 
emission function is sufficiently homogeneous to contribute to the correlation 
function. Thus HBT measures "regions of homogeneity" in the source and 
their variation with the momentum of the particle pairs. As we will see, the 
latter is the key to their physical interpretation. 



3.2 YKP parametrization for the correlator and HBT radius parameters 

A full characterization of the source in terms of its second order space-time 
variances requires knowledge of the 10 parameters (x^x^). These quantities 
appear in the expression (|13| ) for the correlation function but this expression 
still uses all four components of the relative momnetum q^. However, as al- 
ready noted only three of the four components are independent, due to the 
mass-shell constraint (ll7j). Thus only 6 linear combinations of the variances 
{x^x v ){K) are actually measurablea 

If the source is azimuthally symmetric around the beam axis, this counting 
changes as follows: Even if the source is azimuthally symmetric in coordinate 
space, the emission function S{x, K) in phase space is for finite K no longer 
azimuthally symmetric because the transverse components K ± of the pair 
momentum distinguish a direction transverse to the beam direction. There 
remains, however, a reflection symmetry with respect to the plane spanned 
by K and the beam axis. If we call the direction orthogonal to this plane y, 
all mixed variances which are linear in y must vanish due to this reflection 
symmetry, and the correlator must be symmetric under q y — > —q y - Thus 
only 7 non-vanishing variances survive in general, of which, due to the 

mass-shell constraint ( |T7| ) only 4 linear combinations are measurable. 
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Before the correlator (]2^) can be fit to experimental data, the redundant 
components of q must first be eliminated from the exponent of the Gaussian 
via (0). We use a cartesian coordinate system with the z-axis along the 
beam direction and the x-axis along K±. Then (3 = (f3±,0,{3i). We assume 
an azimuthally symmetric source (impact parameter ss 0) and eliminate from 
(B8T) q x and q y in terms of q\ — q 2 x + q%, qi and q° . This yields the YKP 
parametrizationE 



C(q, K) = 1 + exp 



-R 



Iql-Rj (qf - (q°f) (i? 2 + flf) (q ■ U) 



(29) 



Here R±, Rn, Rq, U are four JT-dependent parameter functions. U(K) is a 
4- velocity with only a longitudinal spatial component: 

U(K)= 7 (K)(l,0,0,v(K)), with 7 (30) 

V 1 — v z 

Its value depends, of course, on thii measurement frame. The "Yano-Koonin 
velocity" v{K) can be calculated in an arbitrary reference frame from the 
second central space-time moments of S(x, K). It is, to a good approximation, 
the longitudinal velocity of the fl uid e lement from which most of the particles 
with momentum K arc cmittcdEiliJ For sources with boost-invariant longi- 
tudinal expansion velocity the YK-rapidity associated with v(K) is linearly 
related to the pair rapidity 

The other three YKP parameters do not depend on the longitudinal ve- 
locity of the observer. (Thisdistinguishes the YKP form ( ^9|) from the Pratt- 
Bertsch parametrizationlllZlll which results from eliminating q° in (p8|).) Their 
physical interpretation is easiest in terms of coordinates measured in the frame 
where v{K) vanishes. There they are given byl§l 

Rl(K) = (f), (31) 

R\{K) = ((z-(A//3^) 2 )-(/V/^)V) *(S 2 ), (32) 

Rl(K) = ({t-S:/f3 1 _f)~(f)/^^{P). (33) 

R±, i?|| and i?o thus measure, approximately, the (JT-dependent) transverse, 
longitudinal and temporal regions of homogeneity of the source in the local co- 
moving frame of the emitter. The approximation in ( |32"| , [33"|) consists of dropping 
terms which for the model discussed below vanish in the absence^ transverse 
flow and were found to be small even for finite transverse flowu-3 Note that 
it leads to a complete separation of the spatial and temporal aspects of the 
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source. This separation is spoiled by sources with (x 2 ) ^ {y 2 }- For our source 
this happens for non-zero transverse (in particular for large) transverse flow 
r/f, but for opaque sources where particle emission is surface dominated^ this 
occurs even without transverse flowed 



4 A model for a finite expanding source 

For our quantitative studies we used the following model for an expanding 
thermalized sourceEl 



Mxcoshfa-F) 
1 ' ' 8^ 4 At 6P 



K-u(x) (r-r ) 2 r 2 (r) - r] ) 2 



T[x) 2(At) 2 2i? 2 2(A?7) 2 



(34) 



y 2 , the spacetime rapidity r\ — \ ln[(i + z)/(t — z)], and the 



longitudinal proper time r = \Jt 2 — z 2 parametrize the spacetime coordinates 
x 11 , with measure d 4 x — rdr drjr dr d<f). Y = i 1h[(Ek + K£)/ (Ek — Kl)} and 
M± = y/ m 2 + Kj_ parametrize the longitudinal and transverse components of 
the pair momentum K. \/2R is the transverse geometric (Gaussian) radius 
of the source, tq its average freeze-out proper time, At the mean proper time 
duration of particle emission, and A77 parametrizes the finite longitudinal ex- 
tension of the source. T(x) is the freeze-out temperature; if you don't like the 
idea of thermalization in heavy ion collisions, you can think of it as a param- 
eter that describes the random distribution of the particle momenta at each 
space-time point around their average value. The latter is parametrized by a 
collective flow velocity u^{x) in the form 

u^{x) = (cosh i] cosh f]t(r), sinh rj t (r) e r , sinh 77 cosher)) , (35) 

with a boost-invariant longitudinal flow rapidity r\i = r\ [v\ = z/t) and a linear 
transverse flow rapidity profile 

Vt(r) = Vf (^) • (36) 

i]f scales the strength of the transverse flow. The exponent of the Boltzmann 
factor in (|34|) can then be written as 

K ■ u(x) — Mj_ cosh(Y — 77) cosher) — K.ye. r sinh % (r) . (37) 

For vanishing transverse flow [rjf = 0) the source depends only on Mj_, and re- 
mains azimuthally symmetric for all K± . Since in the absence of transverse flow 
the /3-dependent terms in (p2|) and (B3) vanish and the source itself depends 
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only on M±, all three YKP radius parameters then show perfect Mj_-scaling. 
Plotted as functions of M±, they coincide for pion and kaon pairs (see Fig. [y, 
left column). For non-zero transverse flow (right column) this Mj_-scaling is 
broken by two effects: (1) The thermal exponent (pST] ) receives an additional 
contribution proportional to K± = >/ j^T m2 ■ (2) The terms which were ne- 
glected in the second equalities of ( p2| , p3| ) are non-zero, and they also depend 
on (3± = K±/Ek- Both effects induce an explicit rest mass dependence and 
destroy the Mj_-scaling of the YKP size parameters. 

5 It-dependence of YKP parameters and collective flow 

Collective expansion induces correlations between coordinates and momenta 
in the source, and these result in a dependence of the HBT parameters on the 
pair momentum K. At each point in the source the local velocity distribution 
is centered around the average fluid velocity; two points whose fluid elements 
move rapidly relative to each other are thus unlikely to contribute particles with 
small relative momenta. Essentially only such regions in the source contribute 
to the correlation function whose fluid elements move with velocities close to 
the velocity of the observed particle pair. 

5. 1 The Yano-Koonin velocity and longitudinal flow 

Fig. U shows (for pion pairs) the dependence of the YK velocity on the pair 
momentum K. In Fig. ||a we show the YK rapidity Y YK = | ln[(l +«)/(l — «)] 
as a function of the pair rapidity Y (both relative to the CMS) for different 
values of K±, in Fig. ||b the same quantity as a function of K± for different 
Y. Solid lines are without transverse flow, dashed lines are for rjf = 0.6. For 
large K± pairs, the YK rest frame approaches the LCMS (which moves with 
the pair rapidity Y) ; in this limit all pairs are thus emitted from a small region 
in the source which moves with the same longitudinal velocity as the pair. For 
small K± the YK frame is considerably slower than the LCMS; this is due to 
the thermal smearing, of the particle velocities in our source around the local 
fluid velocity it M (:r)li-a The linear relationship between the rapidity Y YK of the 
Yano-Koonin frame and the pion pair rapidity Y is a direct reflection of the 
boost-invariant longitudinal expansion flowli-3 For a non-expanding source Y YK 
would be independent of Y. Additional transverse flow is seen to have nearly no 
effect. The dependence of the YK velocity on the pair rapidity thus measures 
directly the longitudinal expansion of the source and cleanly separates it from 
its transverse dynamics. 

The NA49 data for 160 A GeV Pb+Pb collisions MM show very clearly 
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M ± (MeV/c 2 ) M ± (MeV/c 2 ) 



Figure 1: The YKP radii Rj_, ijy, and Ro (top to bottom) for zero transverse flow (left 
column) and for i)j = 0.6 (right column), as functions of Mj_ for pairs at Y cm = 0. Solid 
(dashed) lines are for pions (kaons). The breaking of the Mj_ -scaling by transverse flow is 
obvious in the right column. For nonzero transverse flow Ro also does not agree exactly with 

the effective source lifetime \J (t 2 ) (lower right panel). Source parameters: T = 140 MeV, 
Arj = 1.2, R = 3 fm, r = 3 fm/c, At = 1 fm/c. 
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800 K ± (MeV/c) 



Figure 2: (a) The Yano-Koonin rapidity for pion pairs, as a function of the pair c.m. rapidity 
Y, for various values of K± and two values for the transverse flow rjf. (b) The same, but 
plotted against for various values of Y and r]f. Source parameters as in Fig. Ill 



such a more or less linear rise of the Yano-Koonin source rapidity with the ra- 
pidity of the pion pair Th is confirms, in the most transparent way imaginable, 
their earlier conclusionl^PJ based on the Y-dependence of the longitudinal radius 
parameter Ri in the Pratt-Bertsch parametrization that the source created in 
200 A GeV S+A collisions expands longitudinally in a nearly boost-invariant 
way. 

It should be noted that this longitudinal flow need not be of hydrodynam- 
ical (pressure generated) nature. In a description of nuclear collisions as a 
set of longitudinally oriented and independently fragmenting nucleon-nucleon 
strings, the string fragmentation process would also lead to a strong correla- 
tion between the longitudinal positions and momenta of the created hadrons. 
Thus a similar linear rise of the YK-rapidity with the pair rapidity would be 
expected in jet fragmentation (with the z-axis oriented along the jet axis). It 
would be interesting to confirm this prediction in e + e~ or pp collisions. 



5.2 M±- dependence of YKP radii; transverse flow 

If the source expands rapidly and features large velocity gradients, the "regions 
of homogeneity" contributing to the correlation function will be small. Their 
size will be inversely related to the velocity gradients, scaled by a "thermal 
smearing factor" y/T/M± which characterizes the width of the Boltzmann 

distributionB If one evaluates the expectation values ( |jT[]33| ) by saddle point 
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integration one finds for pairs with Y = 



R \ 


— R* : 


Rq 


= (At*) 2 




= Ll, 



with 



111 

m ~ r 2 + r 



(38) 
(39) 
(40) 



1^ + — > (41) 



flow 

2 



(At*) = (AtY+2^/^+LI-to) , (42) 
1 11 



Ll (r A V ) 2 L% 



(43) 



ow 



where Rn w and Lfl ow are the transverse and longitudinal "dynamical lengths 
of homogeneity" due to the expansion velocity gradients: 



^ m±) = ^fe=%(W^' (44) 



L ^ M ^ = T ^wr^ik> (45) 

where ui is the longitudinal 4-velocity. 

Thus, for expanding sources, the HBT radius parameters are generically 
decreasing functions of the tran sverse pair mass M± . The slope of this decrease 
grows with the expansion ratelU'O (this cannot be seen in the saddle point 
approximated expressions above). Longitudinal expansion affects mostly the 
longitudinal radius parameter i?y and the temporal parameter i?oP-3 the latter 
is a secondary effect since particles from different points are usually emitted 
at different times, and a decreasing longitudinal homogeneity length thus also 
leads to a reduced effective duration of particle emission (see lower panels in 
Fig. |l|). The transverse radius parameter R± is invariant under longitudinal 
boosts and thus not affected at all by longitudinal expansion (upper left panel 
in Fig. |l|). It begins to drop as a function of Mj_, however, if the source 
expands in the transverse directions (upper right panel) . Comparing the lower 
two left and right panels in Fig. Hi one sees that the sensitivity of i?|| and Rq to 
transverse flow is much weakerlU Transverse (longitudinal) flow thus affects 
mostly the transverse (longitudinal) regions of homogeneity. 
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While longitudinal "flow" is not necessarily a signature for nuclear col- 
lectivity but could be "faked" as discussed at the end of the previous sub- 
section, transverse flow is much more generic in this respect: there is clearly 
no transverse collective dynamics in the ingoing channel, and the only mecha- 
nism imaginable for the creation of transverse flow is multiple (re-)scattering 
among the participants and secondaries, leading ultimately to hydrodynamic 
transverse expansion. 

Unfortunately, the observation of an Mj_-dependence of R± by itself is 
not sufficient to prove the existence of radial transverse flow. It can also be 
created hv^ther types of transverse gradients, e.g. a transverse temperature 
gradient To exclude such a possibility one must check the Mj_-scaling 
of the YKP radii, i.e. the independence of the functions Ri(M±) (i =_L, ||,0) 
of the particle rest mass (which is not broken by temperature gradients). Since 
different particle species are affected differently by resonance decays, such a 
check further requires the elimination of resonance effects. 

6 Resonance decays 

Resonance decays contribute additional pions at low M± ; these pions originate 
from a larger region than the direct ones, due to resonance propagation before 
decay. They cause an M j -dep endent modification of the HBT radii. 

Quantitative studies BjHl2 have shown that the resonances can be subdi- 
vided into three classes with different characteristic effects on the correlator: 

(i) Short-lived resonances with lifetimes up to a few fm/c do not propagate 
far outside the region of thermal emission and thus affect R± only marginally. 
They contribute to Ro and R\\ up to about 1 fm via their lifetime; R\\ is larger 
if pion emission occurs later because for approximately boost-invariant expan- 
sion the longitudinal velocity gradient decreases as a function of time. 

(ii) Long-lived resonances with lifetimes of more than several hundred fm/c 
do not contribute to the measurable correlation and thus only reduce the cor- 
relation strength (the intercept at q = 0), without changing the shape of the 
correlator^ The reason is that they propagate very far before decaying, thus 
simulating a very large source which contributes to the correlation signal only 
for unmeasurably small relative momenta. 

(iii) There is only one resonance which does not fall in either of these two 
classes and can thus distort the form of the correlation function: the u> with 



"A reduced correlation strength in the two-particle sector could also arise from partial 
phase coherence in the sourceEd By comparing two- and three-particle correlations, the 
intercept reducing effects of resonances can be eliminated, fw*d the degree of coherence resp. 
chaoticity in the source can be unambiguously determinedCJ 
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Figure 3: The inverted q-variance R± and the kurtosis Ax (the index s in the figure stands 
for "sideward") at Y = as functions of K±. Left: r)j = (no transverse flow). Right: 
r)j = 0.3. The difference between dashed and solid lines is entirely dominated by ui decays. 
Source parameters as in Fig. hi except for R = 5 fm. 



its lifetime of 23.4 fm/c. It contributes a second bump at small q to the cor- 
relator, giving it a non-Gaussian shape and complicatingliil the extraction of 
HBT radii by a Gaussian fit to the correlation function. At small M± up to 
10% of the pions can come from u> decays, and this fraction doubles effectively 
in the correlator since the other pion can be a direct one; thus the effect is not 
always negligible. 

In a detailed model studylll we showed that resonance contributions can 
be identified through the non-Gaussian features in the correlator induced by 
the tails in the emission function resulting from resonance decays. To this end 
one computes the second and fourth order g-moments of the correlatorl!^ The 
second order moments define the HBT radii, while the kurtosis (the normalized 
fourth order moments) provide a lowest order measure for the deviations from 
a Gaussian shape. We found that, at least for the model (|34|), a positive 
kurtosis can always be associated with resonance decay contributions (Fig. |[ 
left panel). Strong flow also generates a non-zero, but small and apparently 
always negative kurtosis (Fig. ||, right panel). Any Mj_-dependence of R± 
which is associated with a positive M^-dependent kurtosis must therefore be 
regarded with suspicion; an Mj_-dependence of R± with a vanishing or negative 
kurtosis, however, cannot be blamed on resonance decays. 

In our model, the first situation is realized for a source without transverse 
expansion (left panel of Fig. 1): At small M± the u> contribution increases 
R± by up to 0.5 fm while for M± > 600 MeV it dies out. The effect on R± 
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is small because the heavy u) moves slowly and doesn't travel very far before 
decaying. The resonance contribution is clearly visible in the positive kurtosis 
(lower curve) . For non-zero transverse flow (right panel) there is no resonance 
contribution to R± ; this is because for finite flow the effective source size for the 
heavier u> is smaller than for the direct pions, and the w-decay pions thus always 
remain buried under the much more abundant direct ones. Correspondingly 
the kurtosis essentially vanishes; in fact, it is slightly negative, due to the weak 
non-Gaussian features induced by the transverse flow. 



7 Opaque sources 

The emission function ( |34| ) is only one of an infinity of possible source parame- 
trizations. It is chosen in such a way that it allows easy implementation of most 
features which we believe are important for the sources created in heavy ion 
collisions. There is, however, one important physical situation which cannot 
be parametrized in any reasonable way by the formula (|34|) : if the source emits 
particles not from the entire volume, but only from a thin surface layer. This 
is how the sun radiates photons, and this is also an often suggested picture for 
the case that a QGP is created in the collision which then hadronizes slowly in 
a deflagration-type strong first order transition by surface emission of hadrons 
from the edge of the QGP blob. 

The significance of such a phe nomenon for HBT interferometry was real- 
ized by Heiselberg and VischerE3 who pointed out that an effective emission 
region which is part of a thin surface layer has a much smaller extension in 
the "outward" or x-direction than in the "sideward" or y-direction. In other 
words, such "opaque sources" have (x 2 — y 2 ) < 0. Depending on the degree 
of opacity (the thickness of the surface layer relative to the source_radius) this 
difference can be large and negative. The authors pointed out that this 
leads to the possibility of a smaller "outward" than "sideward" HBT radius 
parameter in the Pratt-Bertsch parametrization, even at K± = 0. Recently 
B. Tomasik showed \=-li that in the YKP parametrization opacity effects would 
show up even more spectacularly by a "lifetime parameter" R 2 , which would 
diverge to — oo in the limit Kj_ — > resp. (3± — ► (see Eq. (p3|)). 

The source (|4|) can be made opaque by multiplying it by the factor! 

ex p(- K l -f) 
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Figure 4: The YKP parameter Rq(M±) for midrapidity pairs for various combinations of A 
and rjf. Solid (dashed) lines are for pions (kaons). Thin lines: A = 1 fm = R/7; medium 
lines: A = 7 fm = _R; thick lines: A = 14 fm = 2R. Left: rjj = 0; right: rjf = 0.6. Source 
parameters: T = 100 MeV, Ar] = 1.3, t = 7.8 fm/c, A T = 2 fm/c, R = 7 fm. 



where A is the mean free path, 

, . _ v 2 f°° _,' 2 , 

feflt = 'eff(?", (p) = e / e dx with y = r sin0, x = r cosc/> . 

J X 

(47) 

is the effective travelling distance of the emitted particle through matter in 
the source (|34|), and k = y/8/% is a parameter which is adjusted in such a way 
that particles emitted from the cejiter suffer the same suppression as in the 
model of Heiselberg and VischerE3 who use a box distribution instead of the 
Gaussians in (|34|). 

Fig. ^ shows the interesting YKP parameter, the "temporal" radius pa- 
rameter Rq, as a function of M± for sources with different degrees of opacity. 
With or without transverse flow, the crucial features of opacity are clearly vis- 
ible: the negative contribution ~ (x 2 — y 2 ) in (|33|) drives Rq to negative values 
at small K±, and this happens the sooner the shorter the mean free path A, 
i.e. the thinner the surface layer is. 

Comparing the solid curves for pion pairs in Fig. ^ with the Pb+Pb dataE^ 
in Fig. H below I conclude that mean free path values A < R are essentially 
excluded. (This conclusion has in the meantime been checked by more exten- 
sive parameter studies.) In other words, the sources created in 160 A GeV 
Pb+Pb collisions are not at all opaque, but rather "transparent" , meaning 
that particles are emitted from the whole volume by bulk rather than surface 
dominated freeze-out. 
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8 Analysis of Pb+Pb data 



In Fig. ^| we show a numerical fit of the YKP radius parameters, using the 
expressions (|3l|)-(|33|) with our model source j3^), to data collected by the 
NA49 collaboration in 158 A GeV/c Pb+Pb collisionsr^r^l Please note that 
this fit refers to only a single rapidity slice of the available data, and it does 
not include resonance decays (although we do not expect the latter to change 
things much, except for reducing tq by about 1 fm/c, the average lifetime of 
the shortlived resonances, see Sec. ^). The fit result must therefore be taken 
with some care. A comprehensive simultaneous analysis of all single particle 
spectra and two-particle correlation data from Pb+Pb collisions is in progress. 

After At7=L3 has been adjusted to reproduce the width of the pion rapid- 
ity distributional the parameters To and Arj are essentially fixed by the magni- 
tude of and Ro . The radius R is fixed by the magnitude of R± (K± = 0) once 
the temperature T and transverse flow rjf are known. The M^-dependence of 
R± fixes T and rjf, albeit not independently: essentially only the combination 
rjf y/M±/T, i.e. the velocity gradient divided by the thermal smearing factor, 
can be extracted^^l This is similar to the single particle spectra whose M±_- 
slopes determine only an effective blushifted temperature^ T c g = T^J j^j^- 
The correlations between T and rjf are, however, exactly opposite in the two 
cases: for a fixed spectral slope T must be decreased if rjt increases while a 
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fixed Af^-slope of R± requires decreasing values of rjf if T is reducedlfHI The 



combination of single-particle spectra and two-particle correlation thus allows 
for a separate determination of T and rjf. 

For the fit in Fig. [|the freeze-out temperature was set by hand to T = 100 
MeV. The resulting flow parameter ijf— 0.6 corresponds to an average trans- 
verse flow velocity v/=0.58. This combination of T and rjf results in single- 
particle spectra with roughly the right shape. Somewhat higher temperature 
values of around T — 120 MeV as advocated by Kampferf^j with an average 
transverse flow velocity Vf = 0.43 c produce a somewhat flatter decrease of R± 
with Mj_ but still appear to be consistent with the data inside the error bars. 
Fitting R± with even higher temperatures results in larger ^/-values which 
leads to single particle spectra which are much too flat. 

Let us discuss in more detail the numbers resulting from this fit. First, the 
transverse size parameter R=7 fm is surprisingly large. Resonance contribu- 



te 



The transverse flow 



tions are not expected to reduce it by more than 0.5 fmL 
correction to R± is appreciable, resulting in a visible transverse homogeneity 
length of only about 5.5 fm at small K±, but even this number is large. R=7 fm 
corresponds to an r.m.s. radius r rms = y/ (x 2 + y 2 ) w 10 fm of the pion source, 
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Data: Pb+Pb 158 AGeV, 
NA49 preliminary, 
Y ab = 4-4.5, FLCMS, 



Calculation: 



rj f = 0.6, R = 7.0 fm, 

x = 7.8 fm/c, Ax = 2 fm/c, 

Arj = 1 .3, T = 1 00 MeV 
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Figure 5: R±, Ru and R for 158 A GeV/c Pb+Pb collisions 



pair momentum. The dataare from the NA49 Collaboration^ 



s functions of the transverse 
The lines are a fit with the 



model (|3J), with fit parameters as given in the figure. 
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to be compared with an r.m.s. radius rf^ s = 1.2x A 1 / 3 * a/2/5 fm = 4.5 fm for 
the density distribution of the original Pb nucleus projected on the transverse 
plane.]] This implies a transverse expansion of the reaction zone by a linear 
factor 10/4.5 = 2.2. That we also find a large transverse flow velocity renders 
the picture consistent. The longitudinal size of the collision region at the point 
where the pressure in the system began to drive the transverse expansion can 
be estimated as follows: for the source to expand in, say the y-direction from 
VVy. n . t . al = 1.2AV3/V5 fm = 3.2 fm to vW fina i = R = 7 fin with an 
average transverse flow velocity of at most v = 0.58 c (the freeze-out value 
determined from the fit of R±) requires at time of at least (7 — 3.2)/0.58 fm/c 
= 6.5 fm/c. Due to the selfsimilarity of the longitudinal expansion the longi- 
tudinal dimension of the source grows linearly with r. If the total expansion 
time until freeze-out is given by the fit parameter tq — 7.8 fm/c, the source 
expanded in the 6.5 fm/c during which there was transverse expansion by a fac- 
tor 7.8/(7.8-6.5) = 7.8/1.3 — 6 in the longitudinal direction. We conclude that 
the fireball volume must have expanded by a factor 6 * 2.2 2 w 30 between the 
onset of transverse expansion and freeze-out! This is the clearest evidence for 
strong collective dynamical behaviour in ultra-relativistic heavy-ion collisions 
so far. 

The local comoving energy density at freeze-out can be estimated from the 
fitted values for T and r\f. The thermal energy density of a hadron resonance 
gas at T = 100 MeV and moderate baryon chemical potential is of the order of 
50 MeV/fm 3 . The large average transverse flow velocity of (vr) w 0.58 implies 
that about 50% flow energy must be added in the lab frame. This results 
in an estimate of about 0.050 GeV/fm 3 x 1.5 x 30 w 2.2 GeV/fm 3 for the 
energy density of the reaction zone at the onset of transverse expansion. This 
is well above the critical energy density 6 ^ < 0.9 GeV/fm 3 predicted by lattice 
QCD for deconfined quark-gluon mattertj Whether this energy density was 
fully thermalized is, of course, a different question. It must, however, have 
been accompanied by transverse pressure (i.e. some degree of equilibration 
of momenta must have occurred already before this point), because otherwise 
transverse expansion could not have been initiated. 

9 Conclusions 

I hope to have shown that 

'in my lecture given at the wotjkshop I made an embarrassing error by a factor \/2 which 
is also contained in the writeupEj of my talk at Hirschegg in January 1996. I would like to 
thank D. Ferenc for pointing out this error. 
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• two-particle correlation functions from heavy-ion collisions provide valu- 
able information both on the geometry and the dynamical state of the 
reaction zone at freeze-out; 

• a comprehensive and simultaneous analysis of single-particle spectra and 
two-particle correlations, with the help of models which provide a real- 
istic parametrization of the emission function, allows for an essentially 
complete reconstruction of the final state of the reaction zone, which can 
serve as a reliable basis for theoretical back-extrapolations towards the 
interesting hot and dense early stages of the collision; 

• simple and conservative estimates, based on the crucial new informa- 
tion from HBT measurements on the large transverse size of the source 
at freeze-out and using only energy conservation, lead to the conclusion 
that in Pb+Pb collisions at CERN, before the onset of transverse ex- 
pansion, the energy density exceeded comfortably the critical value for 
the formation of a color deconfined state of quarks and gluons. There 
is, however, no evidence for long time delays due to hadronization of the 
QGP, and pion freeze-out occurs in bulk rather than from the surface 
of the collision fireball. This is in line with lattice results which pre- 
dict at most a weakly first order confinement transition, and with other 
evidencelH for rapid hadronization. 
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